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Resonances in 28 Si+ 28 Si. II 

Analyses for the Angular Distributions and Angular Correlations — 
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Resonances observed in the 28 Si + 28 Si collision are studied by the molecular model. In 
the preceding paper, it is clarified that at high spins in 28 Si + 28 Si (oblate- oblate system), the 
stable dinuclear configuration of the system is equator-equator touching one, and that the 
axially asymmetric shape of the stable configuration of 28 Si + 28 Si gives rise to a wobbling 
motion (if- mixing). There, the normal modes around the equilibrium have been solved 
and various excited states have been obtained. Those states are expected to be the origin 
of a large number of resonances observed. Hence their physical quantities are analyzed 
£N| . theoretically. The results are compared with the recent experiment performed in Strasbourg 

£Nj ■ and turn out to be in good agreement with the data. Disalignments between the orbital 

angular momentum and the spins of the constituent 28 Si nuclei in the resonance state are 
clarified. Moreover the analyses of the angular correlations indicate characteristic features 
for each normal-mode excitation. Thus it is possible to identify the modes, and a systematic 
experimental study of angular correlation measurements is desired. 



Subject Index: 223 

§1. Introduction 



Intermediate resonances observed in heavy-ion scattering have offered intrigu- 



ing subjects in nuclear physics] i n th e preceding paper,® (hereafter referred as 
paper I), the authors have studied dinuclear molecular structure of the 28 Si+ 28 Si sys- 
tem. In the present paper, we investigate physical quantities by using the molecular 
I/"") \ wave functions and compare the results with the experiments. 

Betts et al. firstly observed a series of resonance-like enhancements at 6 cm = 90° 
in elastic scattering of 28 Si + 28 Si, in the energy range from E\ a \j = lOlMeV to 
128MeV with broad bumps of about 2MeV width. They gave spin assignments of 
J = 34 — 42 by the Legendre-fits to the elastic angular distributions for each bump, 
which correspond to the grazing partial waves.®'® They further closely investi- 
gated angle- averaged excitation functions for the elastic and inelastic scatterings, 
and found, in each bump, several sharp peaks correlating among the elastic and 
inelastic channels.®'® The total widths of those resonances are about 150keV, and 
the inelastic decay strengths are enhanced and stronger than the elastic one, which 
suggests that they are special eigenstates of the compound system. Similar sharp 
resonance peaks are observed by Zurmiihle et al. in the 24 Mg + 24 Mg system.® 
In those systems, the decay widths of the elastic and inelastic channels up to high 
spin members of the 24 Mg or 28 Si ground rotational band exhaust about 30% of the 
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total widths, whereas those into a-transfer channels are much smallerPM These 
enhancements of symmetric-mass decays strongly suggest dinuclear molecular con- 
figurations for the resonance states. It is also noted that the widths of the elastic 
channel are rather small, for example, a few keV, which is quite different from high 
spin resonances in lighter systems such as 12 C + 12 C and 16 + 16 0. 

Taking into account the difference from resonances in lighter systems and the 
level density of the sharp resonances observed, we have developed a new dinucleus- 
molecular model for the high spin resonances in the 24 Mg + 24 Mg and 28 Si + 28 Si 
systems 113.113 JED >D3} in which two incident ions are supposed to form a united com- 
posite system. It rotates as a whole in space with internal degrees of freedom orig- 
inating from interaction of the deformed constituent ions. This is in contrast with 
the viewpoint of the band crossing modelP^ 1 

In the paper I™ we have applied the 
model also to 28 Si + 28 Si, though there 
is a distinct difference between prolate- 
prolate and oblate-oblate systems. In 
the former, the stable configuration is 
pole-to-pole one due to the prolate de- 
formation of 24 Mg, while in the latter, it 
is the equator-to-equator configuration 
due to the oblate deformation of 28 Si. 
Therefore, the former composite system 
is axially symmetric in the equilibrium, 
while the latter is triaxial. Then, in 
the latter, strong X-mixing is kinemati- 
cally induced, and results in a wobbling 
motion Thus we extend the molecular 
model to the axially-asymmetric config- 
urations. In practice, we do not treat 
Coriolis terms in the hamiltonian ex- 
plicitly, but we diagonalize the hamil- 
tonian of the asymmetric rotator due to 
the triaxiality to obtain new low-lying 
states. They are expected to correspond 
to the sharp resonance peaks within 
each bump of the grazing J. There- 
fore, it appears that the present model 
is promising for the sharp high spin res- 
onances observed in 28 Si + 28 Si, but it is difficult to proceed further, i.e., to assign 
which mode of excitation corresponds to which peak observed, due to less detailed 
experimental information. Fortunately, there are exceptional data on the resonance 
at E cra = 55.8MeV in 28 Si + 28 Si, as we will discuss in detail in the present paper. 

A new facet of the resonance states has been explored by the 28 Si + 28 Si scatter- 
ing experiment on the resonance at E cm = 55.8MeV at IReS StrasbourgP Fi gure 1 
shows those elastic and inelastic angular distributions, where the solid lines in the 
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Fig. 1. Experimental angular distributions for 
the elastic and inelastic scattering for 
28 Si + 28 Si at £ CM = tt.&MeV^ m Solid 
curves show L = 38 Legendre fits for com- 
parison. 
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lower three panels show L = 38 Legendre fits. The oscillating patterns in the elastic 
and inelastic channels 2 + , (2 + , 2 + ) are found to be in good agreement with L = 38, 
which suggests L = J = 38 dominance in the resonance, namely, disalignments of 
the fragment spins and the orbital angular momentum. 7-rays emitted from the 
emerging fragments, the excited 28 Si nuclei, have been also measured with Air detec- 
tors in coincidence with two 28 Si fragments detected at = 90°. Those angular 
correlation data show characteristic "m = 0" patterns in the axis normal to the re- 
action plane, which suggests that the fragment spins are in the scattering plane and 
is consistent with the disalignments observed in the fragment angular distributions. 
The disalignments have never been observed before in heavy-ion reactions, which 
indicates that the 28 Si + 28 Si system is completely different from 12 C + 12 C etc. and 
even from the 24 Mg + 24 Mg system which exhibit spin alignments.^'^®'^'^ It, 
thus, is worth to emphasize here that the angular correlation data provide crucial 
information on configurations of the constituent nuclei in the resonance state, which 
could not be obtained otherwise. Therefore, we concentrate our analyses on the 
data on the resonance E cm = 55.8MeV in the following. With the new data on the 
angular correlations, we are now able to select "which mode is really a candidate for 
the resonance". 

For this purpose, we have developed a method of analysis for " particle-particle-7 
angular correlations", in the framework of R- matrix theoryP^'^l 1 The results ob- 
tained by the analysis of the angular correlations show that X-mixing due to the 
rotational motions of the triaxially-deformed system is found to be indispensably 
important for understanding the "m = 0" dominance of the angular correlations, 
which brings novel aspects of "wobbling motion"^ in the resonances, or in the reac- 
tion mechanism. 

In §2, formulation of the i?-matrix theory is briefly reminded, and expressions of 
angular correlations are given in the framework. In §3, theoretical analyses and com- 
parisons with the experimental data are made. In §4, we discuss the mechanism of 
the spin disalignments and the observation conditions of the resonances. Conclusions 
are given in §5. 

§2. Formulation on decay properties and angular correlations 

For calculations of the decay properties of the molecular resonance states, we 
prepare necessary expressions, based on .R-matrix theory,^ since we have the energy 
spectra and the wave functions of the molecular states obtained in the paper IP 1 

2.1. R-matrix formalism 

We write the channel spin wave functions ip as 

^t(i 1 i 2 )im i = Yl {hhMiM 2 \IM I )xr l hM 1 ■ W2M2, (2-1) 

K'hM 2 

where (I, Mi) denote the channel spin, and XnhMi describe the states of the sep- 
arated constituent nuclei. With (Tj/jMj), Tj specifies the internal state of the i-th 
constituent nucleus, and {hM,i) denote the spin quantum numbers. Since we con- 
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sider excitations to the members of the ground-state band of 28 Si nucleus, Tj is taken 
to be the same as that of the ground state, and thus we omit it in the following 
descriptions. Complete channel wave functions are given as 

%,lM z im = Ud J- C Yi m (9,</))lf)( Ill2 ) JMl , (2-2) 

\/VrX r 



where u c i{r c ) denotes the radial wave function between two constituent nuclei in the 
channel c, which specifies I2), I being the orbital angular momentum. For the 
analyses of resonance states with the total spin (J, M), hereafter we use generalized 
channel wave functions in the representation {(Iityll, JM} as 

y c iuM = 2 (IlMMJM^j^YUOA). (2-3) 

Mi,m 

For the radial wave function u c i(r c ), the asymptotic functions are given with 

u i~) =e i* cl{Gd _ iFd)? 

u {+) = e - i ^(G d + iF cl ), (2-4) 

where u , ' and denote incoming and outgoing waves, respectively, with a c \ 
being Coulomb phase shift. F c i and G c \ denote the regular and irregular Coulomb 
wave functions. 

Specifying the incident channel with two 28 Si nuclei in the ground state, with 
I = J and m = with z-axis parallel to the beam direction, the wave function of the 
system for the external region is written as 

& ~ A- y c iuM - > J Uc'i'i'cii—pL — y C 'i'vjM, (2-5) 

where c denotes the initial channel (Ii = I2 = 0), and d denotes the final channel, 
respectively. By using i?-matrix formula,^ we obtain the collision matrix, 

TT — U cl \k c ,ac) x \ ^ u \c'I'l'^XcIl f f-\ 

U c >ri',cii - -j+y- -6c' c -«2^ N J(E-W J Y { ' 

u d '(k c ,a c ) A jv a<^ rV\) 

where a c denotes channel radius. The second term of Eq. (|2-6p is a sum of the 
contributions from the resonance states A. W x is a resonance pole, i.e., W x = 

-A - V A ■ ;l "< 1 -\\ ' '■' ' ' 

l \cll 



— f-Tjf , and N-[ corresponds to a factor for the normalization of the resonance 
state A, which is close to unit. u{ n are defined by 



j y/2k c a c j , 

u \di = -rry- -Txclb (2-7) 

u y cl '{k c ,a c ) 

with the reduced widths TacJZ f rom the amplitudes of the resonance states in the 
channel c, and Jxdl are §i ven by 
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where denotes decaying resonance states, S being the surface at the channel 

radius a c . Note that notations u J XcIl correspond to the definition of Eq. (|2-7|) . but 

with the incoming waves uL instead of u^ 1 " , as is usual in the i2-matrix theory. 



Calculations of Eq. (|2-8p need relations between molecular wave functions and the 
generalized channel wave functions, which are described in Appendix A. 

For inelastic scattering, U^pi' c n is given only by the second term. We apply 
one level approximation, i.e., we consider the energy region close to a resonance 
level A and replace the sum of the second term of Eq. (|2-6p by one term with the 
resonance state A, the effect of the neglected terms being taking into account by the 
modification of the hard sphere scattering phase shift as is explained later, and then 
we have 

TT — „■ U \c'I'l' U XcIl (o n\ 

it, Nf(E-Wf) V ) 
We obtain cross sections for the reaction d 7^ c, 

I q . rj I 2 I . tj ^ I 2 

_ _ _ n sr^ \ TT 12 _ . v \ u xc'i'i'\ \ a \cii\ f0 in \ 

9 ^h l M - SJ ^k\^n(E-E(f +i nr <240) 

where gj = (2 J + l)/(2Ji + 1)(2J2 + 1) is the statistical factor. Now we define 
partial widths in the channel c' of the level A by 



I'V I'V 
where P c i denote so-called penetration factor given by 



and the total width of the level A is obtained by the sum of the partial widths, 

A J = E r Ac- (2-13) 

c 

Hence, with \N^\ ~ 1, we obtain the Breit-Wigner one- level formula, 

vr Ac' Ac ro -tA\ 

"""•"mB-W + iif}- ( ' 

In particular, for "on resonance", i.e., for E = E{, we obtain 



u<„„Mm = (2.i5) 

X X 

where the indications of J and A for the resonance state are omitted. 4> c [ denote hard- 
sphere scattering phase shifts, which are defined by tan</> c j = F c i(k c , a c )/G c i(k c , a c ) . 
Note that in the expression Eq. f|2- 13|) for the total width, the channels c should be 
taken generally not only for channels such as the elastic, the single and mutual 2 + but 
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also for those associated with higher fragment spins. Furthermore experimentally 
it is suggested that the contributions from symmetrical-dinucleus channels exhaust 
only about one third of the total widths.® Therefore we practically adopt the total 
width r( = 150keV for Eq. (f2T5j) which is suggested by the experiments P 

To calculate "{{ell by fl2-8p . we adopt the model wave functions obtained 
by the bound-state approximation for . Relations between the molecular wave 
functions and the channel wave functions necessary for the calculations are given 
in Appendix A. In order to avoid large dependence on the channel radius due to 
gaussian damping of the model radial wave function, which is due to the harmonic 
approximation, we have corrected tails of the radial wave functions by smoothly 
connecting G c {s in the concerning channels. At E cm = 55.8MeV the distances of 
the connecting points are in the range of R c = 8.3 — 8.8fm, and each channel radius 
a c is taken to be larger than the connecting distance, which means pure Coulomb 
interaction is assumed outside a c . From R c to a c , internal waves in the channel 
(c, Z) of are assumed to behave to be proportional to (G c i + iF c i), like the 

outgoing waves. Hence the channel radius dependence is removed except for the 
small renormalization effect due to N^, which is defined by = J W r ^ M *W^ dr 
with the integration inside a c . By using wave functions normalized in the bound 
state approximation, is usually real and close to unit. In detail, a small difference 
occurs from modified radial tails around a c , which, of course, depends on E cm and 
a c . The difference in each channel appears to be a few percent at the maximum, 
and for simplicity we have approximated the value of N\ to be equal to one in the 
present calculations. 

For the elastic scattering in which various partial waves contribute from I = to 
larger than grazing partial ones, it is better to rewrite Eqs. (|2-5p and (|2-6p explicitly 
with I\ = 1% = and Z's. We have 



-L-j2(2i + i)i l 

^/vkr ^ 



e m Fl + _( e 2m_ Ul)u M 



Pi(cos6), (2-16) 



with c being dropped in the notations, and the differential cross section is given by 

2 

V(2Z 4-lVf 2 '"' - TL)Pi(c,nsf)\ 

2k 



da 

1q 



/a(«) + iE( 2l + 1 )( e2w, - p ') P '( cosfl ) 



i 



(2-17) 



where fc(9) denotes the Coulomb scattering amplitude, and Ui denotes the collision 
matrix for the Z-partial waves. As for the collision matrix, we modify Eq. (|2-6p as 
follows. Firstly in order to do one-level approximation, we divide the second term of 
Eq. (12 -6p . a sum of the contributions from the resonance states A, into two groups, 

U t = "£ (Mc) - i V j^St Ja - j _ , (2-18) 



where the third term is the contribution from the resonance under consideration, and 
the second term is those from the other levels. In general, the second term includes 
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not only the contributions from the resonance states with the 28 Si + 28 Si configu- 
rations, but also those from the resonance states with the different configurations 
and from the compound nucleus states. Since the contributions to the total widths 
from the symmetrical dinuclear channels exhaust only about one third, as already 
mentioned, the contributions from the other channels are significant. The effects of 
the second term brought from the other decay channels and the fusion channels in 
average are expected to appear as absorption in the elastic channel. Thus we replace 
the sum of the first and second terms of the r.h.s. of Eq. (|2- 18|) with nuclear scat- 
tering amplitudes which is phenomenologically proposed by the study of heavy-ion 
scattering,^ and obtain the collision matrix for the /-partial waves as 

Ui = Ai-i "ft*** j. , (2-19) 

where A[ denotes the elastic scattering amplitude for the Z-partial wave. Supposing 
that the lower partial waves than the grazing one with l gr ~ J would be strongly 
absorbed as usual in the heavy-ion scattering, we assume reflection coefficients to 
be gradually increasing as I comes to be over the grazing one, i.e., we apply smooth 
cutoff modelf® 

e 2i(cr ; +<5 ; ) 

M = l + exp[-(Z-(J + 4))/A]' (2 ' 20) 
with nuclear phase-shifts being b\ = 25 /{l + exp((Z — (J + A))/A{)}. We put A = 2, 
A[ = 3.3 and 5 = 0.1 to reproduce the experimental data in the region of 9 = 20°— 60° 
where the differential cross section is very rapidly decreasing as 9 increases 

Finally it is noted that for the collisions of identical particles, scattering cross 
sections are modified due to the symmetry between the particles. For the final chan- 
nels with identical particles (the cases of I\ = I2), the expression of the differential 
cross sections is given by \(l/y2){f c >(9, 4>) + f c '(^ — 9, <ft + 7r )}\ 2 X 2, where f c > denote 
the amplitudes of the outgoing waves on the spherical surface. The factor (I/1/2) is 
due to the normalization of the incident channel and the multiplied factor 2 corre- 
sponds to the detections of the scattered particle and the recoil particle. Thus the 
scattered partial waves are allowed only for I = even, and the cross sections a c i c 
in Eq. (|2T4p are modified as well, by factor 2. Furthermore careful calculations of 
d(j{9)/dQ by using Eqs. (|2T7p and (|2-19p are necessary. The resonance term is ob- 
tained with the wave functions with boson symmetry. On the other hand the other 
amplitudes fc(&) e tc. are not symmetrized. Thus the differential cross section for 
the elastic scattering of the spinless identical particles is given as 



da 

7q 



1 I 11^ V ^* 

{f (6) + /(vr - 6)} -^{2J + l) ArJ/ ^ \ ir j ) PA<xx 1 



T2 V W ^ nn " )} " 2~k^ J ^ l) NliE- W J 



x 2, (2-21) 



with the amplitude f(9) by the Coulomb and nuclear scattering except for the single 
resonance term, 

f(9) = fc(9) + ^ 5> Z + l)(e 2MTi - A,)^(cob 9). (2-22) 
l 
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2.2. Fragment-fragment-^ angular correlation 

For the calculations of the angular correlations, we need the scattering ampli- 
tudes in magnetic substates (mi,ra2) representation, which are connected to the 
subsequent 7-ray emission process. The amplitudes for the mutual 2 + excitation are 
given by 

X mim2 (k',k) = V (22m 1 m 2 \I'm 1 + m 2 ){I'l'm 1 + m 2 m'\JM) 

IK Z — ' 

I'l'm' 

x Ul, v Y] M (k)Y Vm ,{k'), (2-23) 

where k and k' denote the initial and final relative momenta between two 28 Si nuclei, 
respectively, and the incident (elastic) channel assignment for the collision matrix is 
omitted. Here we have assumed a single resonance with a total angular momentum 
J. For the single excitation, of course, we have a similar expression as the above, 
by putting I 2 = and m 2 = into the CG coefficients of the mutual-channel spin 
coupling in Eq. (|2-23|) . The transition amplitudes for the 7-ray emissions from the 
polarized nuclei are discussed by several authors For two photon emissions from 
the mutual excitation, the amplitudes are proportional to the scattering amplitudes 
and the strength of the photon emission as 

A hh~ Xrn^ik^kXOOlH^lhm^OOlHazlhrnz), (2-24) 

mim,2 

where (00\H a \Iim,i) denote the transition matrices for 7-ray emissions, which give 
transition rates, i.e., by well-known perturbation theory, 2-ir/h- KOOl-ffo-liimi)! 2 /}^ 
for the direction dfl of the emission, with p being the level density of the final states. 
In the case of the present angular correlations, photons are rapidly emitted with the 
half life of 700fs, which means almost all the first 2 + states of 28 Si nuclei finish their 
transitions into the ground states just after the decays of the resonance compounds. 
Thus the problem which we will discuss is not the magnitudes of the transition 
rates but angular distributions of the 7-ray intensities over 4-7r-detectors, which were 
measured by Eurogam Phase II. The 7-ray intensity distributions are given by 

(00|# a (Im) I Im) ~ D 1 ^ (0 7 , 7 , 0) (JJm - m|00) (0| \T(I) | \I) 

(_-i\I-m 

~ ^=^dL(# 7 )e- m< \ (2-25) 

where a denote right/left-hand circular polarizations of the emitted 7-rays, i.e., 
a = ±1. After sum over them for the square of the absolute values of Eq. (|2-24[) . we 
obtain 7-ray angular correlations, 



Whi a (0i,<h;62,<h;k',k) = £ |^ lff2 ' 2 



Viia 

<71<T2 

o\ (72 mi "12 m'i Tn'n 



x (OOlff^IimiXOOli^JIim'O* 

x (00\H a2 \I 2 m 2 )(00\H a2 \I 2 m' 2 )* , (2-26) 
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where fa) denote the directions of the emitted photons. In the experiment, only 
one of the two emitted photons is detected in most cases, even with the EurogamP^ 
Therefore we take up the distribution of one photon after averaging over the distri- 
bution of the other photon. By taking average over the angles of the second photon 
(#2, fa), the last line of Eq. (|2-26p is reduced to 8 m2m i due to the integration for the 

phase factor e _ ^ m2_m 2^ 2 appearing from Eq. (|2-25p . Hence the Air 7-ray intensity 
distribution is given by 

w hh (0i,fa;k',k)^ {T / x ^m 2 (k',k)X* m[m2 (k',k)} 

mi raj m 2 

x (OOlfl^l/imiJCOOlIT^lJimi)*, (2-27) 

where sum over 777,2 is associated with the average on one photon. 

Now with the expression of Eq. (|2-27p . we are able to compare the theoretical 
results with the experimental Air 7-ray distribution. As is usually done, it is useful 
to determine the contributions from the magnetic substates.^ 1 In order to obtain 
the probabilities in the magnetic substates of the fragments, we integrate Eq. (|2-27j) 
over fa (averaged around the z-axis). Then only diagonal elements with mi = 777^ 
appear because of the phase e - J ( m i- m i)^i from Eq. (|2-25j) . Thus the normalized 
angular correlations are expressed as follows: 

w(e 1 ) = Y,P m w m {e 1 ), (2-28) 

m 



E 
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Fig. 2. 7-ray intensity distributions W m (8-y) 
of the pure magnetic substates for m = 
0, ±1 and ±2. From the top panel, m — 0. 
\m\ = 1 and \m\ = 2 are shown, respec- 
tively. 



(a) 
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Fig. 3. 7-ray intensities for particle-particle- 
7 angular correlations with the maximum- 
aligned configuration J — 38 with I' = 
34(1' = 4) of the mutual 2+ channel. For 
the quantization axes in three panels, see 
text. 
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with 

P m ~^|X mm2 (fc / ,fc)| 2 , (2-29) 

and 

Wm(fcr) = J E {V5/^d 2 ma (9^)} 2 , (2-30) 

where P m denote probabilities in the magnetic substates m, and W m (9-y) denote E2 
7-ray angular distributions of m-substates. The intensity distributions, W m (S 7 ) are 
well-knownJ23 i.e., W (9) = (15/8vr) sin 2 9 cos 2 0, W ±1 (6) = (5/16vr)(l -3cos 2 # + 
4cos 4 9) and W±2{9) = (5/167r)(l— cos 4 9). In Fig. 2, they are shown for convenience, 
from the top panel Wo(9), W±\{9) and W±2(9), respectively. 

IReS 4ir 7-data were given for the probabilities in the magnetic substates, where 
three different quantization axes are taken: (a) the beam axis, (b) the axis normal 
to the scattering plane, and (c) the axis perpendicular to (a) and (b) axesP^ 1 Note 
that (b) is important for understanding spin alignments, where the z-axis is normal 
to the scattering plane. Later in Fig. 8, the "m=0" pattern is seen in (b), which 
suggests fragment spins are in the scattering plane. In contrast, in Fig. 3, we see 
the "m=2" pattern clearly in (b), where theoretical results of the maximum-aligned 
configuration (J = 38 with I' = 34(1' = 4)) are displayed, in three panels (a)~(c) 
corresponding to three quantization axes. 

We specify the initial beam direction k in the coordinates corresponding to 
each quantization axis. The final fragment directions k' are determined by the 
experimental procedures p3>C3 [ n which two large-area position-sensitive particle 
detectors are located symmetrically on either side of the beam axis in the horizontal 
plane, to take triple coincidence of fragment-fragment-7. The 28 Si fragments are 
detected in the angular range of 9 cra = 90° ± 7° in the reaction plane, and with 
the vertical acceptance ±4°. Thus the final fragment direction is approximately 
corresponding to (c)-axis. 

For quantitative analyses of the angular correlations, it is important to take the 
angular range of particle detectors into account. One may think that the size of the 
particle detector is enough small to neglect the effect. But in extremely high spin 
resonances such as L ~ 38, amplitudes of the partial waves rapidly oscillate along 9 cm 
in the reaction plane; it needs only 2° from the maximum to zero amplitude as is seen 
in Fig. 1. Hence we need to take the average over the detector area, especially along 
9 cia . Note that such an effect is indispensably important for the angular correlations 
in the single excitation channel, because due to the Bohr condition we have no 
m = ±1 component at exact 9 cm = 90° p^ 1 while actually we have rather large 
m = ±1 components as seen in (a)-axis of the experimental datcP^ (see Fig. 11). 
On the other hand, the vertical acceptance ±4° is not serious. Since the orbital 
angular momentum vector is almost perpendicular to the reaction plane, the waves 
very slowly vary in the vertical direction, though they propagate rapidly oscillating 
in the reaction plane. Therefore we need to integrate over only one dimension, i.e., 
along # cm with interval 90° ±7°. To be sure, we have compared the results between 
those with one-dimensional integration and with two-dimensional one. Actually the 
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differences are very small, less than 1% of the total amounts. 

To compare theoretical results with experiment, background effects in the angu- 
lar correlation measurements should be considered.^ Since the excitation yields in 
the single and mutual 2 + channels are seen to have backgrounds of over 50% some 
contributions are expected, probably from non-resonant aligned configurations. As 
is shown in Fig. 3, the aligned configurations appear to provide \m\ = 2 contributions 
in (b), while they provide dominant \m\ = 1 contributions together with significant 
m = ones in (a) and (c). Note that angular correlations with the aligned configu- 
ration for the single excitation appear to be almost the same as those for the mutual 
excitation shown in Fig. 3. 

2.3. Internal wave functions: molecular model of the 28 Si + 28 Si system 

Internal wave functions for the ZZ-matrix calculations are, of course, those of 
dinuclear molecular model, the characteristic features of which are briefly reminded 
here, since they are given in detail in the paper IK« 

Interaction between two nuclei is described with internal collective variables, 
i.e., the orientations of the poles of the constituent nuclei in the rotating molecular 
frame. Assuming a constant deformation and the axial symmetry of the constituent 
nuclei, we are dealing with seven degrees of freedom, 

(qd = (0 1 ,e 2 ,0 3 ,R,a,p 1 ,fo), (2-31) 

as illustrated in Fig. 4, where (i?, #2j#i) is the relative vector of the two 28 Si. As 
internal degrees of freedom, the orientations of the symmetry axes of the two 28 Si are 
described with the Euler angles (aj,/3j) which refer to the molecular axes. a.\ and 
a2 are combined into #3 = (a\ + a2)/2 and a = {a\ — a^jl. All those orientation 
dependences of the interaction are described by a folding-model potential. 

It is found that at high spins the dinuclear system with oblate-deformed con- 
stituent nuclei has the equilibrium in an equator-equator touching configuration with 
the parallel principal axes, and that the relative distance between the two 28 Si nuclei 
is 7 — 8fm indicating a nuclear compound system with hyperdeformation. The barrier 




Fig. 4. The dinuclear configuration and the coordinates in the rotating molecular frame for an 
oblate-oblate system. The molecular z'-axis and the seven degrees of freedom of the system 
are displayed, where the ol\- and a2-degrees are to be combined into 6*3 = (ai + 012) /2 and the 
degree of twisting a — (ai — a2)/2. The figure is taken from Ref.[3|. It is the same as published 
in Refs.E} and 03]). 
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position (or saddle point) is 9 — lOfm, greatly outside from that of the usual optical 
potential. Molecular ground state configurations are well bound by the barrier up 
to J = 40. This theoretical maximum spin is in accord with the bumps observed in 
the grazing angular momenta. 

Since the interaction potential does not couple the states of different if-values, 
at first, we assume the eigenstates of the system to be the rotation-vibration type 
with a good i<C-quantum number, 

Vx~D J MK (0i) X K{R,a,fa,fa), ( 2 - 32 ) 

then the problem to be solved is of internal motions associated with the variables 
(R, a, fa, fa). Couplings among various molecular configurations are taken into ac- 
count by the method of normal mode around the equilibrium configuration, which 
gives rise to the molecular modes of excitation. We expand the effective potential at 
the energy minimum point (the equilibrium relative distance R e = 7.6fm, a = tt/2 
and fa = fa = vr/2), and adopt the harmonic approximation to obtain the normal 
modes; the results are as follows. The radial motion has no coupling with the other 
angle variables, and it is an independent mode in itself. The motions associated with 
the /3-degrees are well confined to be vibrational ones, which are classified into new 
modes, butterfly: fa = (A fa + Afa)/y/2 and anti-butterfly: fa = (A fa - Afa)/\[2 
around a = tt/2, respectively. As for the a-degree, the confinement in the present 
folding potential appears to be unexpectedly weak, and hence the motion is close to 
a free rotation, which we call twisting rotational mode. Thus we write the internal 
wave functions as 

XK(R,a,fa,fa) = f n (R)4> K (a)<pt + (f3+,o>)ip~_(fa,a), (2-33) 

where ip^ + (fa, a) and ip~_ (/?_ , a) denote the functions describing the fa and fa 
modes, respectively, a's in tp+ ,(/3+,a) and ip~_(fa,a) are due to the strengths of 
the confinements in the /3±-motions, which vary together with the a-motion. 
The eigenenergy of the system is given as 



E J (n,n + ,n_,K,(u,7r a )) = E (R c ) + 

+ {n + ^)^R 



h 2 \J{J + 1)- K 2 -1 + K 2 



fiR 2 21 



+ (n + + n_ + l)^ + £"(7T a ), (2-34) 

where the energy is specified by the quantum numbers (n, n+, n_, K, (z/, ir a )), with v 
as a dominant frequency of the a-motion and with the parity 7r a about the reflection 
with respect to a = tt/2. The first and second terms on the r.h.s. of Eq. (|2-34p 
are constant energies that are given by the interaction potential and the centrifugal 
energy at the equilibrium, respectively, where I denotes the moment of inertia of the 
constituent nuclei, 28 Si, the value of which is estimated from the excitation energy 
of the 2f state. The term (n+ + n_ + l)fouo gives vibrational energies for the fa 
motions without the a-dependence, and finally E"(7r a ) is the energy for the a-motion 
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including the a-dependent contribution to the energy of the /3-motions. The values 
of the vibrational energy quanta for the butterfly and the anti-butterfly modes are 
both about 4MeV, but the excitation energies of those modes appear twice, 8MeV, 
since states of K = even with one vibrational quantum are not allowed due to 
the boson symmetry. The excitation energy is close to that of the radial excitation. 
Although the excited state of the radial mode is not bound in the present calculations, 
the possibility of the radial-mode resonance is not completely excluded, because it is 
likely that the interaction between two 28 Si would be more attractive than the present 
folding potential with the frozen density approximation. It is noted that the in- 
dependence in <fn± (/^±) a ) is rather tedious in calculations of the partial widths in the 
next section. Hence we dropped the a-dependence in them, i.e., to be <p n ±(l3±) with 
the average value 4MeV of huip, which means that the eigenenergies of the butterfly 
and anti-butterfly modes are taken to be degenerate. However the properties of those 
modes are essentially retained with no problem. Of course, for exact calculations, we 
need total wave functions with the parity and boson symmetries, which are described 
in Appendix B of the paper some examples of the wave functions of the normal 
modes being also given there in Appendix C. 

Wobbling motion (K -mixed states) 

One of the characteristic features of 
the spectrum obtained theoretically is a 
series of low-energy if-rotational exci- 
tation due to axial asymmetry around 
molecular z-axis, which is in contrast 
with the 24 Mg + 24 Mg casePDffl) One 
can understand the reason immedi- 
ately from Fig. 5, where the upper 
conflguration( 24 Mg + 24 Mg) has axial 
symmetry as a total system, but the 
lower one for 28 Si + 28 Si has axial asym- 
metry. 

A triaxial system preferentially ro- 
tates around the axis with the largest 
moment of inertia. By the definition of 
the axes in the lower panel of Fig. 5, 
we have the moments of inertia as Ix > 
Iy » Iz, due to the nuclear shape. Thus the system, which is seen as two pancake- 
like objects ( 28 Si's) touching side-by-side, rotates around X-axis normal to the reac- 
tion plane. Such a motion is called wobbling, and K is not a good quantum number. 
Namely, we expect that the eigenstates are if-mixed. 

Regarding the 28 Si + 28 Si system as a triaxial rotator, as it is usual for poly- 
atomic molecules, we diagonalize the hamiltonian with an inertia tensor of the axial 
asymmetry, which gives rise to mixings of if -projections of the total spin JpD The 
resultant motion should be called as "wobbling mode"P2) The energy spectrum is 
displayed in Fig. 6(b), compared with the spectrum without if -mixing in Fig. 6(a). 




Y 



Fig. 5. Equilibrium configurations of two di- 
nuclear systems. The upper portion is 
for 24 Mg — 24 Mg and the lower one for 
28 Si- 28 Si. The figure is taken from Ref.0). 
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Fig. 6. Energy spectra of the 28 Si + 28 Si system for J = 38. (a) Molecular normal modes without 
/("-mixing, (b) After /("-mixing, with indications of the modes under the levels. The figure is 
taken from Ref. 



Now the states of low lying iT-series are not the eigenstates by themselves, but are 
recomposed into new states. It is much interesting that we again obtain several 
states including the K = component as a result of K-mixing, which should show 
up themselves in the scattering. Those states are closely located in energy and so in 
good agreement with several fine peaks observed in the experiment. 

As an analytical prescription, in the high spin limit (K/J ~ 0), the solution is 
a gaussian, or a gaussian multiplied by an Hermite polynomial, 



F n {K) 



exp 



1 (K 
'2\T 



2i 



The width parameter b is given by 

b=(2J 2 I K /A) 1 /\ 



(2-35) 



(2-36) 



which depends on the ratio between the non-axiality (the difference between Ix and 
Iy) and Iz as 



1 

A 



2\ly 



1 

Ix 



1 



1 



- — — 

2Uy + Ix 



(2-37) 



To calculate angular correlations we use those analytic forms in Eq. (|2-35|) , which 
is simple and intuitive way to understand the extent of if-mixing. Of course we can 
utilize numerical values obtained in the diagonalization procedure, but the values 
are almost the same as those given by the analytic form. For the lowest state Fq{K) 
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of Eq. ()2-35p . we have the wave function for the wobbling ground state as 

H M -J2 ex ^- K2 / 2b ^ D MK(^)XK(R,a,p 1 ,(3 2 ), (2-38) 

K 

where it should be emphasized that, in general, xk can be any molecular mode of tri- 
axial deformations, such as the ground-state configuration (parallel equator-equator 
one), the butterfly mode and the anti-butterfly mode. Then the spin orientations of 
two 28 Si nuclei are expected to be in the plane, consistent with "m = 0", because the 
nuclei rotate around the axes perpendicular to their symmetry ones. The magnitude 
of b estimated by Eq. (|2-36p is 1-85, for example, for the values of the moments of 
inertia used in the calculations for the energy spectrum in Fig. 6. This is the largest 
value expected, because we assumed a static configuration there, in which the zero- 
point motions of the twisting and butterfly modes are neglected. Thus, a value of 
b = 1.5 is adopted for the calculations of the physical quantities in the next section. 
There, it is also examined to what extent a change of b affects the spin disalignments 
with the orbital angular momentum. 

As for the relation between the molecular model and the asymmetric rotator 
intuitively introduced, we mention that it is clarified in §4.2 of the paper by 
using a simple example of the dinuclear system of "one deformed nucleus and a 
spherical nucleus" , and thus we do not repeat here. 

§3. Results 

3.1. Structures and partial decay widths of the resonance states 

In order to analyze structures of the resonance states by the molecular model, we 
have estimated their partial widths. In Table I, the partial widths of the molecular 
normal modes with spin J = 38 are given, up to the mutual 2 + channel. The results 
for the 4 + channels are not shown, since the experimental angular distributions of 
those channels do not exhibit resonance behaviors fES-HZ) Magnitudes of the widths 
are calculated with the theoretical level energies, and are obtained to be in the range 
of several keV to several tens keV, which is consistent with experiment. For corn- 
Table I. Theoretical estimates for partial widths (in keV) in the elastic, single 2 + and mutual 2 + 
channels of the molecular resonances with spin J = 38 for the theoretical level energies, which 
are obtained in §2.3. For the wobbling strength, b — 1.5 is adopted. In parenthesis, widths are 
given with the assumption of the resonance energy 55.8MeV. 



Molecular states 


Energy 


r ei 






molecular ground state 


51.5MeV 


2.1 


1.2 


0.16 






(19) 


(13) 


(2.2) 


twisting 


54.4MeV 


0.36 


0.24 


6.2 






(0.69) 


(0.46) 


(13) 


butterfly 


58.9MeV 


28 


1.5 


21 






(8.3) 


(0.43) 


(5.0) 


anti-butterfly 


60.3MeV 


40 


1.5 


69 






(8.3) 


(0.25) 


(15) 
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parison, the magnitudes of the widths are also calculated with the assumption that 
the energy of each state is shifted to the observed resonance energy of 55.8MeV, and 
are given in parenthesis. Among the elastic and inelastic channels, characteristic 
features are seen for each normal mode. As for the molecular ground state, the 
elastic width is larger than those in the excited channels, which is inconsistent with 
the experimental characteristics of the narrow resonances. Those weak excitations 
of the theoretical prediction are inferred to be due to the weak confinement in the 
a-degree of freedom, which permits almost free rotation for the solution with the 
folding potential. In the butterfly and twisting modes, the mutual 2 + channel shows 
strong excitation, but the single 2 + excitation is very weak, and thus the charac- 
teristics are inconsistent with experiment. Therefore neither of the normal modes 
is completely consistent with the characteristics of the enhanced excitations seen at 
E cm = 55.8MeV. 

It is, however, meaningful to investigate with a stronger confinement in the 
a-degree of freedom for the molecular ground state, because, in the touching con- 
figuration, the confinement possibly increases due to induced deformations of the 
constituent nuclei. So we investigate the cases with stronger confinements for the 
a-degree with the harmonic oscillator, by using fkj a as a parameter. This means 
that we introduce a stronger confinement of the gaussian distribution around the 
equilibrium, instead of the internal-rotation-like motions, although how much the 
strength of the confinement should be increased cannot be predicted by the folding 
model within the frozen density approximation, but is to be quantitatively investi- 
gated with the polarization effects or by dynamical treatments of deformations of 
the constituent ions. The polarization effects of the ions expected to appear will be 
discussed later in §4.3. 

Table II. Theoretical estimates for partial widths (in keV) in the elastic, single 2 + and mutual 
2 + channels for the molecular ground state with spin J = 38 at the experimental resonance 
energy E cm — 55.8MeV. Rather confined configurations for the molecular ground state are 
inspected for (case A: hw a = 2MeV, ftwp = 4MeV), (case B: hw a — ftwp = 4MeV) and (case C: 
hui a = hwp = 12MeV). For the wobbling strength, b = 1.5 is adopted. 



States 


r ei 




■* (2+,2+) 


molecular ground state A 


17 


12 


2.9 


molecular ground state B 


12 


9.2 


4.3 


molecular ground state C 


2.6 


3.5 


3.0 



In Table II, a couple of configurations with stronger confinements are introduced, 
and the estimations of the partial widths are given at the energy 55.8MeV. For 
the configurations named A, B and C, strengths of the confinements are set to 
be hu> a = 2MeV, 4MeV and 12MeV, respectively, where in the last case ftwp = 
12MeV is taken consistently with the large value fru) a = 12MeV. As the confinement 
increases, the magnitude of the partial widths in the elastic channel decreases, while 
excitations to the single and mutual channels are relatively enhanced. With the 
configuration C, characteristics of the partial widths are in good agreement with 
the experimental data.-- Moreover the value of the elastic partial width 2.6keV is 
enough small to reproduce the experimentally suggested value, "a few keV'P Thus 
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the molecular ground state with a well confined configuration is a candidate for the 
narrow resonance at 55.8MeV, which is corroborated by the following analyses on 
fragment angular distributions and fragment-fragment-7 angular correlations. 

3.2. Fragment angular distributions for the resonance at E cm = 55.8MeV 

In Fig. 7, fragment angular distributions are displayed for the elastic scatter- 
ing, the single 2 + and the mutual 2 + excitations. Theoretical angular distributions 
with the configuration C of the wobbling ground state are compared with the recent 
datap^* where their magnitudes are normalized to the single excitation. For the ex- 
citations, constant backgrounds of 40% are assumed, because the background yields 
of about 50% are seen in the experimental data observed by Betts et alP In the 
elastic scattering, by the phenomenological use of the potential scattering term as 
in Eq. (|2-19p . calculations are made with strong absorption of lower partial waves 
than the grazing one to reproduce rapid decrease in the differential cross sections 
toward # cm = 90° JDJD Such a rapid decrease is usually seen in the heavy-ion scatter- 
ing |25},[32]i an( j .j. WQU ^ j-^ Qne Q £ necessar y conditions for the observation of the quite 



5 




50 60 70 80 90 100 110 
9 cm (degree) 

Fig. 7. Fragment angular distributions for the elastic and inelastic scattering for 28 Si + 28 Si at 
-Ecm = 55.8MeV. From the bottom those for the elastic, single 2 + and mutual 2 + are shown. 
Solid curves are for theoretical results of the molecular ground state with wobbling motion 
(configuration C), compared with the experimental dataP^ 
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narrow resonances in the elastic scattering. We see good fits for the oscillations, the 
period of which is close to L = 38 Legendre polynomials. Magnitudes of the cross 
sections are qualitatively in good agreement with the data, although the strengths 
of the excitations are still slightly weak. 

3.3. Angular correlations 

In the decay process from the molecular resonance state, 28 Si- fragments emit 
7-rays due to the transition from the first excited state to the ground state. In 
the experiment, particle detectors are set to catch the fragments in the direction 
perpendicular to the incident beam, and the coincident 7 rays are detected by the 
An system of Eurogam Phase II. Figure 8 displays the angular correlation datap^ 
i.e., 7-ray intensity distributions of the £'2-transition observed in the mutual 2 + 
channel and x 2_ fits shown by dashed lines. Three different quantization axes are 
taken in panels (a) — (c), respectively: (a) beam direction, (b) z-axis normal to the 
scattering plane, and (c) z-axis perpendicular to those of (a) and (b). Since 28 Si 
fragments are detected in the angular range of 6 cm = 90° ± 7° in the reaction plane, 
the z-axis of (c) corresponds approximately to final fragment directions. 



(a) 

1 1 




1 1 1 1 1 


(c) 

1 1 
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Fig. 8. 7-ray intensities for particle-particle-7 angular correlations in the mutual 2 channel. Ex- 
perimental data are given by dots, and x 2_n ts curves are displayed by dashed lines. The theo- 
retical results obtained with the molecular ground state (configuration C: ftw a — hwp — 12MeV) 
are shown by solid lines. For the quantization axes in three panels, see text. 
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Table III. The experimental and theoretical probabilities P m in the magnetic substates m, for 
the mutual 2 + channel of the resonance state with spin J = 38 at _E cm = 55.8MeV. The 
experimental values of P m are obtained by x 2-n ts of the dataP^ Theoretical results are given 
for the molecular ground state, where the first line in (a) shows those for the weak confinement 
(solution with the folding potential), while the second line in (a) does those for the strong 
confinement (configuration C: = hoj a = 12MeV), and in (b) and (c) as well. For the 

wobbling strength, b = 1.5 is adopted. 



Quantization 


Pm=Q 




-Pm = ±l 




-Pm = ±2 




axis 


exp. 


theor. 


exp. 


theor. 


exp. 


theor. 


(a) 


0.30 ± 0.08 


0.51 


0.16 ±0.04 


0.07 


0.18 ±0.05 


0.17 






0.20 




0.06 




0.34 


(b) 


0.46 ±0.05 


0.70 





0.01 


0.27 ±0.02 


0.14 






0.54 




0.18 




0.05 


(c) 


0.14 ±0.05 


0.05 


0.17 ±0.03 


0.08 


0.26 ± 0.04 


0.40 






0.14 




0.17 




0.26 



Theoretical results for the molecular ground state (configuration C of a strong 
confinement) are displayed by solid lines. In panel (b), we see a typical "m=0" 
pattern, which means that the spin directions are in the reaction plane. A weak 
"m=0" pattern is seen also in panel (a), while in panel (c), a slight swelling around 
the center suggests weak "m=2". 

Probabilities P m in each magnetic substate m are listed in Table III. Theoretical 
results of the molecular ground state are given and compared with the experimen- 
tal datap^ where the upper line in each axis case is for a weak confinement case 
(solution with the folding potential), while the second line is for a strong confine- 
ment (configuration C), respectively. In (b), theoretical results reproduce dominant 
"m=0" characteristics, for both configurations of the weak and strong confinements. 
Thus, the "m=0" characteristics of the angular correlations do not depend essen- 
tially upon the strength of the confinement in the a-degree of freedom. \m\ = 2 
contributions in (b) appear slightly weak, compared with the experimental data, but 
this suggests \m\ = 2 contribution from the background.^ 1 As for the values of P m 
in (a) and (c), they are also in good agreements with the data, although they vary 
depending on the strengths of the confinements. 

To clarify properties of the normal mode excitations in the angular correlations, 
we make analyses for the twisting and butterfly modes, which are of interest as 
they exhibit strong excitation to the mutual 2 + channel. In Fig. 9, the angular 
correlations of the twisting mode are displayed with dashed lines, and comparison 
with the molecular ground state (solid lines) is given. Apparently, the twisting 
excitation does not fit the experimental data. The characteristic feature of twisting 
is "m=2" dominance in (a), which is due to the excited rotational motion in the 
a-degree of freedom. 

In Fig. 10, results are displayed with dashed lines for the butterfly mode with 
if-mixing. The results for the anti-butterfly mode appear to be quite similar to 
those of the butterfly mode (not shown here). Again, solid lines are those for the 
molecular ground state for comparison. The butterfly excitations show their own 
characteristics with dominant "m=l" patterns in the panels (a) and (b). In addition, 
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they exhibit "m=2" dominance in (c), due to the butterfly vibrational motions. 
Thus all those theoretical results for the normal mode excitations are much different 
from the data at the resonance energy 55.8MeV, and moreover they show quite 
distinguishing characteristics. 

Note that the configurations of the molecular ground state in Figs. 8, 9 and 
Fig. 10 are not completely the same due to changes of the confinement parameter 
value, and one can confirm differences between the results with different strengths 
of confinements in the a-degree of freedom. In Fig. 9, the most weak confinement is 
taken from the solution with the folding potential, and in Fig. 10 a slightly stronger 
confinement is taken (the configuration A: hu a = 2MeV, hup = 4MeV). All those 
configurations give "m=0" characteristics, and the differences between them are not 
significant. 
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Fig. 9. 7-ray intensities for particle-particle- 
7 angular correlations in the mutual 2 + 
channel. Experimental data are given by 
dots. Dashed lines show theoretical re- 
sults of the twisting rotational mode, with 
solid lines for the molecular ground state 
for comparison. For the quantization axes 
in three panels, see text. 
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Fig. 10. 7-ray intensities for particle-particle- 
7 angular correlations in the mutual 2 + 
channel. Experimental data are given by 
dots. Dashed lines show the theoretical re- 
sults of the butterfly state with if-mixing, 
with solid lines for the molecular ground 
state for comparison. For the quantization 
axes in three panels, see text. 
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The magnitude of wobbling (K- 
mixing) is taken to be b = 1.5, which 
is consistent with the asymmetry of the 
28 Si + 28 Si molecular ground-state con- 
figuration, as it is estimated to be b = 
1.85 in the limit of the strong confine- 
ment (without quantum fluctuations). 
The parameter for the X-mixing is ex- 
amined in the range of b = 1.2 — 2.0 
with the configuration C. The value of 
P m= o in (b) varies from 0.4 to 0.8 with 
the increasing value of b, but for the con- 
figurations with the strong confinements 
the resultant curves are seen to be essen- 
tially the same in those values of b. 

Next we move to the single 2 + ex- 
citation. In Fig. 11, the experimental 
data are displayed, together with x 2 ~ 
fits given by dashed lines. 101 The the- 
oretical results for the molecular ground 
state (configuration C of a strong con- 
finement) are shown by solid lines. Sur- 
prisingly, the experimental data are seen 
to be quite similar to those of the mu- 
tual 2 + channel, and again "m=0" char- 
acteristics in (b) are firstly of interest. 
The present model reproduces the char- 
acteristics well. Theoretical results ex- 
hibit totally good fits with the data, al- 
though the fits in panel (a) are not very 
good, with rather large m = contri- 
bution. However, contributions of the 
background from aligned configurations 
are expected to provide \m\ = 1 compc 
deviations of the fits are not serious. 



(a) 
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Fig. 11. 7-ray intensities for particle-particle- 
7 angular correlations in the single 2 + 
channel. Experimental data are given by 
dots, and x 2 -fits by dashed lines, respec- 
tively. Solid lines show the theoretical re- 
sults obtained with the molecular ground 
state with wobbling motion (configuration 
C). For the quantization axes in three pan- 
els, see text. 

ts in (a) and (c), and thus those slight 



§4. Discussion 



4.1. Structure of the Si — Si molecule 

In the paper I™ interaction between two 28 Si nuclei is described with the internal 
collective variables, i.e., with the orientations of the poles of the constituent nuclei 
in the rotating molecular frame. For the dinuclear system with oblate-deformed 
constituent nuclei, an equator-equator touching configuration is found to be the 
equilibrium at high spins where the principal axes of the constituent nuclei are par- 
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allel. The relative distance between the two Si nuclei is 7— 8fm indicating a nuclear 
compound system with hyperdeformation. The barrier position is 9 — lOfm, greatly 
outside from that of usual optical potentials. Molecular ground state configurations 
are well stable by the barrier up to J = 40. This theoretical maximum spin is in 
accord with the bumps observed in grazing angular momenta. 

Couplings among various molecular configurations are taken into account by the 
method of normal mode around the equilibrium configuration, which gives rise to 
the molecular modes of excitation, such as the radial vibration, the butterfly motion, 
the anti-butterfly motion and so on. The twisting mode {y = 4) appears to be the 
lowest excitation, but the energy may be higher than the present results due to a 
stronger confinement as physically expected, which will be discussed later in §4.3. 

Since the equilibrium configuration has a triaxial shape, we extend our molecular 
model so as to include couplings between states with different X-quantum numbers. 
In practice, we use the asymmetric rotator as an intuitive model, which preferentially 
rotates around the axis of the largest moment of inertia, accompanied with /f-mixing, 
and thus gives rise to the wobbling motion. 

4.2. On the spin disalignments 

An important physical quantity which probes the structure of the resonance 
states is spin alignments of the outgoing particles. The recent experiment on the 
28 Si + 28 Si resonance at E cra = 55.8MeV suggests L = J = 38 dominance in the 
inelastic scattering both to the single and mutual 2 + excitations. Measured particle- 
particle-7 correlations show "m=0" dominance in 7-rays. All these suggest that the 
spin vectors are in the reaction plane 

The analyses are made for the molecular model wave functions. In the state with 
the lowest energy of a given angular momentum J, due to the triaxial configuration, 
the whole system rotates about the axis normal to the reaction plane defined by 
the two pancake-like nuclei. The spins of the 28 Si fragments are thus in this plane, 
since no rotation can occur about the symmetry axes of the 28 Si nuclei. Such a 
property is in agreement with the lack of strong alignments observed in the fragment 
angular distributions. Actually, "m=0" dominances in the angular correlations are 
well reproduced for the single and mutual 2 + excitations. 

Characteristic features of the angular correlations in the normal mode excita- 
tions are found to be dominant "m=2" patterns, i.e., for the twisting rotational mode 
"m=2" appears in (a)-axis, and for the butterfly mode it appears in (c)-axis, respec- 
tively. This means that spin vectors are parallel to the beam direction in the twisting 
mode, while parallel to the fragment direction in the butterfly mode. (The reason for 
those orientations of the spin vectors is explained later.) One may expect that the 
twisting mode and/or the butterfly mode are favorable for the "m=0" characteristic, 
because spin orientations are in the reaction plane. However, more precisely, "m=0" 
means symmetric around the normal axis, which is satisfied by neither of excited 
states such as twisting nor butterfly with a well-defined direction of the spin in the 
reaction plane. The existence of dominant m = 2 substates, together with "m=l" 
patterns in the other panels, are clear characteristics of the normal modes. Thus, it 
is possible to distinguish among the molecular ground state, the twisting excitation 
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and the butterfly excitations, respectively. 

Finally, we discuss the observed spin directions in space, with respect to the 
orientation of the molecular z'-axis at the time of decay, which is related to the 
positions of the particle detectors in the angular correlation measurements. For 
example, in the twisting mode, the 28 Si nuclei rotate around the molecular z'-axis, 
and so the spin vectors are parallel to it. Since the fragments 28 Si are detected 
approximately at cm = 90° (the (c)-axis direction), the molecular z'-axis and the 
spin vectors are expected to be approximately parallel to (c)-axis. However our 
theoretical prediction does not exhibit the "m=2" pattern in (c). Unexpectedly, 
it appears in (a)-axis (the beam axis), as is displayed in Fig. 9 with dashed lines, 
which gives rise to a puzzle. This is an interesting aspect of extremely-high spin 
resonances. Due to the high-speed rotational motion, the final velocity vectors of 
the fragments are considered to be almost perpendicular to the relative vector of 
the fragments in the decay process, like raindrop motions splashed from a rapidly 
rotating umbrella. In order to examine the motions of the 28 Si nuclei after the decay, 
we have analyzed classical orbits and obtained results that for J = 38 the relative 
vector turns round about 70° from the initial angle. As a confirmation, we calculated 
angular correlations with a relatively lower angular momentum J = 8. The results 
show the "m=2" dominance in (c)-axis, as expected by the above explanation. 

4.3. Strengths of confinements 

With respect to the very narrow resonances observed as correlating among the 
elastic and inelastic scattering, partial decay widths are investigated. In the but- 
terfly and the twisting modes, the decay probability amplitudes concentrate to the 
elastic and mutual 2 + channels. This is a characteristic feature of the normal mode 
excitations, which is due to the symmetric motions of the two 28 Si nuclei. Such 
characteristics are, however, not seen in the data at E cm = 55.8MeV. The molecular 
ground state or the radially-excited state obtained by using the folding potential 
shows rather large elastic widths and weak excitations to the inelastic channels. 
Hence, the properties of the partial widths observed at E cm = 55.8MeV have not 
been well reproduced for all those states. 

It is worth to have a close look at the results with the folding-model potential. 
The molecular ground state configuration obtained by using the folding potential 
with the frozen density approximation, exhibits very weak confinement in the twist- 
ing degree of freedom, which results in an almost-free rotation. Intuitively, the ex- 
treme weak or shallow potential in the twisting motion appears unphysical, because 
the two 28 Si nuclei are touching with each other. That is, the folding model with the 
frozen density is not well adequate. Actually, the modification of the confinement 
in terms of larger values of hcu for the twisting and consistently for the butterfly 
motions shows that the partial decay width in the elastic channel greatly reduces, 
while those in the inelastic channels relatively increase. As a result, the magnitudes 
of the partial widths become to be consistent with the experimental characteristics, 
which finally gives good agreements with the experimental data for the fragment 
angular distributions. Thus, the analyses suggest a substantial change of the defor- 
mation of the constituent nuclei and/or formation of a neck between them, which 
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gives rise to a stronger confinement than that obtained by the folding model. This 
is exactly consistent with the existence of the stable triaxial configuration. In ad- 
dition, experimental observations suggesting the stronger confinement are known in 
the 24 Mg + 24 Mg system. Only 30% of the resonance flux appears in the binary 
decay channels of 24 Mg + 24 MgP l A search for the remaining flux to the fusion- 
evaporation channels has been recently made.® Those observations suggest that 
two incident nuclei interact with each other strongly with their density overlapping 
and induced deformations. Their effects would be described effectively by a stronger 
confinement. A mechanism for the strong confinement is intriguing and yet to be 
clarified theoretically. 

4.4. Observation conditions 

Finally, observation conditions of resonances in heavy-ion reactions should be 
discussed. In experimental excitation functions of the 28 Si + 28 Si system which are 
angle- averaged around 8 cm = 90°, background yields of the elastic scattering and 
of the single and mutual 2 + excitations decrease rapidly as the bombarding energy 
increases toward the resonance region.®'^ 1 The same is the case in the 24 Mg + 24 Mg 
system.® Such low background yields are one of the conditions to clearly observe 
resonances. Especially in the elastic scatterings, strong damping of non-resonant 
amplitudes is a necessary condition, because the resonance yields themselves are very 
small. Over those background yields, many prominent peaks would be observed with 
relatively large yields in the channels concerned, which theoretically have to be well 
described by the resonance terms, such as by the i?-matrix theory. Strong absorption 
by the smooth cutoff model has worked as a successful description for the elastic 
background scattering in the angular distributions, as well as appropriate optical 
model analyses. Actually, in the 28 Si + 28 Si system, it is shown that the effects of 
the strong absorption do not necessarily work to eliminate resonance structures, but 
show up very narrow and prominent structures over the very low background. This 
is quite different from the conditions of the resonances observed in lighter systems 
such as 12 C + 12 C and 16 + 16 0, where weak absorption plays a crucial roleP^ 1 

§5. Conclusions 

Resonances in the 28 Si + 28 Si system have been studied by means of the molec- 
ular model, in which the interacting dinuclear system is described by the molecular 
rotation and the internal collective variables for the orientations of the poles of the 
constituent nuclei. The molecular model predicts rotational spectra with a variety 
of intrinsic molecular states, such as the butterfly, twisting and etc. In order to 
explore which one of the various molecular modes agrees with the experimental data 
on the resonance at E cia = 55.8MeV, we have performed comprehensive analyses on 
available physical quantities, i.e., on partial widths, fragment angular distributions 
and fragment-fragment-7 angular correlations. 

Among the various modes, the molecular ground state with J = 38 is successfully 
selected to be a candidate for the resonance observed at E cra = 55.8 MeV. We have 
found that for the resonance state, a nuclear molecule 28 Si - 28 Si of triaxial shape 



Resonances in Si + Si. II 



25 



wobblingly rotates with a very high spin. The present results are the first discovery 
of the modes with JC-mixing in heavy-ion resonance states, though the tilting or 
wobbling mode was once discussed in deep inelastic scattering processes 

In the partial decay widths, each normal mode shows an interesting feature, i.e., 
resonance amplitudes appear to be enhanced in each relevant characteristic chan- 
nels. Moreover, the butterfly modes show fragment spins in the fragment direction 
("m=2"), while the twisting mode does in the beam direction. Since each normal 
mode shows its own distinct characteristics of the angular correlations, it is possible 
to identify each excitation of the modes. Thus, angular correlation measurements 
are a powerful tool for the study of nuclear structures of heavy-ion resonances. Sys- 
tematic measurements are desired, not only on the resonance at E cra = 55.8MeV. 
For example, the same kind of measurements on the other nearby resonances of 
28 Si + 28 Si is strongly called for. 
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Appendix A 

Relation between Molecular Wave Functions and Channel Wave Functions 

for Calculating Decay Properties 

The total system consists of two deformed nuclei, for which we have assumed 
the axial symmetry and constant deformations. We have seven degrees of freedom 
illustrated in Fig. 12(a) except for the center of mass motion for the whole system, 
that is, the relative vector R = (R,02,9i) and the orientations of the deformations 
of the interacting nuclei described with Euler angles (di,/?i) and (012,^2) • Thus 
with those variables, generalized channel wave functions of {(IiI 2 )H, JM} scheme 
are defined as usual in the i?-matrix theory. Specifying the spins I\ and I2 for the 
two constituent nuclei, a function with channel spin / is given by 

^(/i/ 2 )/Mj = ^2 {hhM l M2\IMi)xi 1 M 1 {ai,h)xi2M 2 {pC2,h), (A" 1 ) 
MiM 2 

where XhMiioti, /%) denote internal wave functions for the constituent nuclei 28 Si. 
The generalized channel wave function is given with the orbital angular momentum 
I for the relative motion, 

y(hh)ii,jM = Yl (UMim\JM)S l2 [i> {I 1 I 2 )IM I Ylm(02,0l)], (A- 2) 

Mi,m 
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where S12 = (1 + V12)/ 'y2(l + ^hh) denotes the symmetry operator between two 
28 Si nuclei. The internal wave functions in Eq. (|A-1|) are assumed to be a rotational 
type, such as, 

XliMii&iJi) = \j 2 \^ 2 1 - P A/ i o( a i»^»7»). ( A ' 3 ) 

where -?Q = are given in consistent with the axial symmetry of the constituent 
nuclei assumed, and hence 7, are spurious and do not appear in the l.h.s of Eq. (|A-3p . 

Next we define variables in the molecular coordinate system for the model wave 
functions. As is illustrated in Fig. 12(b), we define a rotating molecular axis z' of the 
whole system with the direction of the relative vector R of the two interacting nuclei. 
In the molecular model, orientations and/or motions of the intrinsic axes of each 
deformed nucleus are referred to the molecular usual. For describing those 

orientations, we introduce new Euler angles (aj., /Si,7i) and (a2,/32,72), relations of 
which with the Euler rotations (a,, ft, 7i) referred to the laboratory system are given 
as 

f2 i {a i: p i ^ i ) = Q , i {a i: p i:li )Q M {0 l ,e 2 ), 1 = 1,2, (A-4) 



(b) Z(e 2 ,e 1 ) 




Fig. 12. The coordinates of an interacting dinuclear system, (a) shows the relative vector R — 
(R,02,9i) and usual Euler angles (on, Pi) of the i-th nucleus referring to the laboratory frame. 
In (b), the molecular z'-axis and the seven degrees of freedom of the system are displayed, where 
the distance R is not indicated explicitly. The third angle 63 is defined by #3 = (ai + 02) /2 to 
give the whole rotation around the z'-axis. The figure is taken from Ref. 
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where and Q[ denote Euler rotations for the i-th constituent nucleus with re- 
spective angles, and Qm denotes rotations of the molecular axes. On the r.h.s. of 
Eq. ()A-4p Q[ denote the successive rotations after Qm', firstly the axes of the i-th 
constituent nucleus rotate up to the directions of the molecular axes by Qm, and 
secondly they rotate referring to the molecular axes by Q[. Since the second rota- 
tions are written as Q[(oti, ft, 7$) = Qm(6i, 6 2 )Qi((Xi, Pi, li)Q~M (61, 62) with rotations 
Qi(ai, /Si,7i) referring to the laboratory system, we obtain 

i7i(5i,/3j,7i) = Q M (0i,0 2 )Qi(ai, $,7*), (A-5) 

and accordingly we are able to define the new Euler angles {oci, ft, 7^) by the relations 

Qi(ai, ft, 7j) = Om (61,92) fki&ijPhli). (A-6) 

Note that 72 are spurious as well as 7$, and so we neglect them in the following 
descriptions. 

The molecular x'-axis would be determined with the two constituent nuclear 
configurations around the z'-axis which are specified by 014. Hence they are combined 
into #3 = (qi + a2)/2 and a = (a± — 02) /2, and then we have molecular coordinates, 

(qi) = (6i,92,9 3 ,R,a,Pi,p 2 ), (A-7) 

where 61,62 and 63 are the Euler angles of the rotating molecular frame with the 
other fours being internal variables. The molecular x'- and y'-axes defined with #3 
are illustrated as x'(9$) and y'(9 3 ) in Fig. 12(b). With the definition of the molecular 
x'-axis with #3, rotational angles of each constituent nucleus should be redefined so 
that 01(6*3) = ai — 6*3 = a and 012(63) = ot2 — 63 = —a. Here after we write them as 
014 simply, and correspondingly the relations Eq. ()A-5p are rewritten as 

Qi(oi,k) = f2M(9i,6 2 ,6 3 )Q i (a i ,p i ). (A-8) 

Now, by means of relations ()A-8p . the internal wave functions of Eq. (|A-3|) are 
written with the molecular coordinates as 

Mi 

We substitute Eq. ()A-9p into Eqs. (jA-ip and (|A-2|) . and use the following relations, 
£ (IihM 1 M2\IM I )D^ i (9 i )D I ^ 2 (9 i ) = (hhWz\I »i + IMfiD 1 ^ 

MxM 2 

(A-10) 

(IlM I m\JM)D I Mit ^ l2 (6 l )Y lm (62,6 1 ) = 

Mi,m 

(II 111 + H2 0| J m + /i 2 )V(2Z + l)/4 7 r J D^ w+M2 (^), (A- 11) 

where 6i is the abbreviation for 61,62,63. Thus we obtain the generalized channel 
wave function expressed with the molecular coordinates, 

3V 2 )«;JM = X)(^0|^)V^^^i2[^(e i )% /3)/Jf (a, J 9i J A)] ) (A-12) 
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with 

= £ (hh^ 2 \IK) e il/a < l0 (/3i)< 2 20 (/3 2 ), (A-13) 

where v = y«x — 1^2 denotes twisting quantum number. As the generalized channel 
wave functions are expanded with the rotational functions for the whole system in a 
series of K = \i\ + fi2, we are able to easily calculate overlapping integrals between 
each channel wave function with {(Iityll, JM} and the model wave functions spec- 
ified with {JMK}. Note that in the i?-matrix formula in §2 and for the generalized 
channel wave functions given above, normalization integrals are defined as usual, 
while the definition of the model wave functions <Pa in Eqs. (|2-32|) and (|2-33p is given 
for the vibrational modes with the volume element dV = dRdadftidfo (see §2.2 of the 
paper for details). Hence, for the internal wave functions in Eq. (|2-8p . we use 
v^D, where W\ denote the model wave functions and D = n^l 2 R 2 I\ sin j3\l2 sin (3 2 , 
square root of which is the additional phase factor of the molecular model wave func- 
tions. 

Finally we mention the property of the channel functions described in the molec- 
ular coordinate system. Due to the angular momentum coupling (IiI 2 iiifJ>2\IK), we 
have restrictions fii < Ji, [12 < I2 and K = [i\ + [12 < I for the summations on K and 
(Ml) M2) i n Eqs. (|A-12p and ()A-13p . For the elastic channel, where I\ and I2 are to be 
equal to 0, the values of [i\ and [12 are zero and so the same as for K and v. Hence 
the elastic channel wave function has no real a-dependence nor ^-dependence, and 
it is described simply by D J MQ {9i,9 2 , #3), i.e., by Yjm(#2> ^1)- Thus, due to the boson 
symmetry, we have only positive-parity states with I = J = even, as usual. 
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